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Abstract 



This paper is concerned with the continuous-time quantum walk 
on Z, Z'^, and infinite homogeneous trees. By using the generating 
function method, we compute the limit of the average probability 
distribution for the general isotropic walk on Z, and for nearest- 
neighbor walks on Z*^ and infinite homogeneous trees. In addition, 
we compute the asymptotic approximation for the probability of 
the return to zero at time t in all these cases. 

1. Introduction 

The concept of quantum walk has its origin in the field of quantum computa- 
tion where the notion of classical random walk has been adapted to the quantum- 
mechanical setting in an attempt to improve the performance of random walk al- 
gorithms. Since its origination in the middle of 90s, this new concept has drawn a 
lot of attention in physical and mathematical literature. 

The early papers that formulated the main ideas of quantum walk are IS and 
|[T4l . From the numerous later papers, we would like to mention ||6l where the 
continuous-time quantum walk was defined and Q which defined and studied the 
discrete-time quantum walk on finite graphs. An introductory review of quantum 
walks can be found in ||9l- For recent developments the reader can also consult 

na. 

In general, a quantum walk is described by a triple {G,ip, Ut), where G is a 
graph, V is a unit-length complex vector function on this graph, i.e., ip T-L = 

(G) (8) C^, IIV'II = 1, and Ut is a family of unitary operators on T-L. 

The interpretation is that the state of a particle at time t is completely described 
by function Utip. Upon measurement at time t, the particle is found at vertex v in 
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State s G {1, 2, . . . , A^} with probability p{v, s, t) = \(Uttp) {v, . (By prob- 
ability we mean a collection of non-negative numbers p{v, s, t) that sum to 1, 

There are two types of quantum walk on graph G. The first type is the discrete- 
time walk. Time is discrete, t S Z, and a step of the quantum walk is given by 
a unitary transformation U, so that Tpt+i = Utl^t- This one-step transformation U 
has some special properties, one of which is locality: Ui^ju = if the graph- 
theoretical distance between vertices u and u is sufficiently large. The discrete- 
time quantum walks on Z and Z*^ have been studied in ifTOl and lU who found that 
their asymptotic behavior is significantly different from the behavior of classical 
random walks. 

In this paper we are going to investigate continuous-time quantum walks (|!6l) 
on infinite graphs. We assume that function ip depends only on the position of the 
particle and time (hence N = 1), that t G M, and that the evolution operators Ut 
are given by the following expression: 

Ut = exp{—iXt), 

where X is a self-adjoint operator ("Hamiltonian") on (G) that respects the 
structure of the graph. One example of such an operator is the discrete Laplacian 
of the graph. We also assume that the initial function i/j is concentrated on one of 
the vertices, the origin of the walk. As usual, the probability to find a particle at 
vertex v if the system is measured at time t is given by \^p {v , t)\'^ . 

The continuous-time walk is not local and the relation between continuous-time 
and discrete-time walks is not yet clear (H). However, the continuous-time has an 
advantage over the discrete-time model in being more tractable analytically. 

For our study, we choose the simplest infinite graphs: the integer lattices Z and 
Z'^, and the homogeneous infinite tree in which every vertex has valency m. 
The continuous-time walk on Z has been previously investigated by several re- 
searchers. In particular, by using the Fourier inversion method, Gottlieb in Q 
derived a formula for the limit probability distribution for a very general class of 
quantum walks on Z. In our paper, an equivalent formula for the limit is obtained 
by using the generating series method instead of the Fourier inversion. The advan- 
tage of this method is that it is easier to generalize it to the case of homogeneous 
trees which is not considered in Q. In addition, this method allows us to show that 
the probability distribution converges to zero faster than any polynomial in time 
in the regions sufficiently fai^ away from the origin. This provides a simple large 
deviation estimate for quantum walks on Z. 
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Let us start with graph G = Z, let the initial function ipi (0) = 60, where 60 
is the delta-function concentrated on 0, and let Ut = exp (iXt) , where X is an 
Hermitian operator on (Z) We are interested in computing 

^Pl (t) = {Si\Ut6o) 

We say that the quantum walk has finite support if there exists a constant L, such 
that Xij = if |i — j\ > L. We call the walk isotropic if Xy = a^-j and a_; = 07, 
where ai are certain constants. If Xij = 6\i_j\_i, we say that the quantum walk is 
nearest-neighbor. 

For a general isotropic quantum walk on Z, define the generating fi^nction of the 
walk by the formula P (9) = o-i^^^^- For example, for the nearest-neighbor 
random walk, P {9) = 2 cos 6. 

Theorem 1.1. Let P (9) be the generating function of an isotropic quantum walk 
with finite support on Z. Assume that ipi (0) = 5q. Let a = l/t and suppose that 
equation P' (9) = —a has K > real solutions 9^ (a) in the interval [0, 27r) . 
Then, the transition amplitude from to I is given by the formula 

(t) = _ 1 y 1 ^uiPie,)+ae,)±m/A + o(-], (1) 



where the sign before iri/A equals the sign of P" (9k) ■ If equation P' {9) = —a 
has no real solutions, then 

for all n and appropriately chosen Cn- 
We will prove this theorem in Section |2l 

Next, define the rescaled probability distribution as p (a, t) = t jV^fat] {t)\ , and 
define the average rescaled probability as 

p{a,T) = ^ [ p{a, t) dt (2) 
V J Jt 

Corollary 1.2. Suppose that equation P' (9) = —a has K > real solutions 
9k (a) in the interval [0, 27r) , and that numbers 

ujk ■■= P{9k) + a9k 

are all different. Then 
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If equation P' (0) = —a has no real solutions, then 

lim p (a, T) = 0. 

T— >oo 

In the case, when some of coincide, formula ^ needs a small adjustment 
which takes into account the positive interference of the exponents with the same 
frequency. 

Proof of Corollary |1.2t From Theorem ll.il it follows that 



t 

After averaging over i, we get 
QED. 

In particular, we obtain the following formula for the probability of return to 
zero. 

Corollary 1.3. Let equation P' (8) = have K > real solutions 9]^ in the inter- 
val [0, 2tt) and assume that P {9k) are all different. Then the average probability 
of the return to zero is 

A somewhat different expression for the limit of the average rescaled probability 
distribution was obtained in Q. Namely, by using the Fourier inversion methods 
Gottlieb obtained the following formula for the limit: 



J{e:-P'{e)ex} 



^d9. 



Here, p^ {X) denotes the limit of the probability to find the particle in the set 
{l]l/t G X} at time t, provided that the initial wave function is ip. The function 
{F*'4)) {9) is defined as Yln&z ''Ane™^- In our case, the random walk is started from 
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the origin and therefore {T*tp) (6) - 
easy to compute 

PjP {[a,a + da]) 



1. Then, by using Gottheb's formula, it is 



Hm ■ 

da 



da 



- E 



1 



IP" I 



-P'(e)=a} 

which is our formula (|3]l. We obtain formula (|3]l by the method of generating 
functions which we will also use for quantum walks on trees. 

For a particular example, consider the nearest-neighbor quantum walk on Z. 
That is, assume that Xjj = Then, P' (9) = — 2sin^ and there are two 

solutions 9k (a) for a < 2: 9i = arcsin(a/2) and 92 = — 9i. Hence, we 
can compute \P" {9k)\ = \/4 — a^, and the average rescaled distribution has the 
following limit: 



p{a,T) 



1 



1 



(4) 



In other words, the average rescaled distribution converges to the arcsine law. (This 
result is the first limit theorem for continuous-time quantum walks, which was 
obtained in ifTTl by using the asymptotics of Bessel functions). 



Figure 1. The limit average probability distribution for the 
quantum walk on Z with P (9) = e"^^^ + e"*^ + e*^ + e*^^. 

Consider another example with P {9) = e~*^^ + e^*^ + e*^ + e*^^. The limit 
probability distribution can be computed numerically by using Theorem 11.11 It 
is shown in Figure [H Note that the points where the probability distribution has 
singularities correspond to local maxima of the function P' (9) = —2(2 sin 20 + 
sin 9). 

Let us summarize the main features of quantum walks on the integer lattice. 
First of all, the length of the interval where the distribution of quantum walk is 
essentially supported is of order t instead of \/t as in the classical case. Then, 
the probability of the return to the origin at time t is of order instead of i^^/^. 
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Finally, unlike the classical case, the limit distribution is not Gaussian and its shape 
depends on the generator of the quantum walk. 

What can be said about the continuous-time quantum walk on Z''? We consider 
here only the nearest-neighbor walk. It turns out that in this case the quantum walk 
on 1,^ factorizes provided it was started from the origin. That is, every transition 
amplitude of the walk on can be written as a product of transition amplitudes of 
the walk on Z. 

For simplicity of notation we consider only the case of 1? . The general case 
is similar. Let (0 be the transition amplitude of the transition from vertex 
(0, 0) to vertex (i, j) . (That is, -(/'(ij) (i) = (%,i)|t^t'^(o,o)) , where (^(o,o) and 
denote the delta-functions concentrated on vertices (0, 0) and (i, j), respectively.) 

Theorem 1.4. Le? the nearest-neighbor quantum walk on 1? be started from the 
origin, i.e., tp (0) = (5(o,o)- Then, 



where ipi (t) is the transition amplitude for the nearest-neighbor quantum walk on 
TL started from the origin. 

We prove this result in |2l Previously, this fact was observed without proof in 
Appendix of HI. 

Let us now turn to continuous-time quantum walks on homogeneous trees. (The 
previous studies of this topic include lH and |[T5l .) We restrict our investigations 
to the case of the nearest-neighbor walk. 

Let G = T„i, the m-valent infinite tree with m > 3, the initial ipbe 6e, where 
e is the root of the tree, and let Ut = exp (iXt) , where X is the adjacency matrix 
of the tree. Let r := 2^/rr^^^. (This is the spectral radius of the operator X.) 
Finally, let us define the following functions of pai^ameter a: 



(*) = V'i {t) ipj {t) 



ui (a) = a arctan 



a 



= + V — 

,2 ^ 



— a 



u)2 (a) = an — uji, 



and 



ifi (a) = — arctan 



m 



a 



vr 



m - 2 ^ 



— a 



,2 



4 



ip2{a) = -iT-ipi. 



Theorem 1.5. Consider the nearest-neighbor quantum walk on the regular infinite 
tree of valency m. Assume ip (0) = 6e, and let ipi (i) be the amplitude of transition 
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from the root to a vertex w which is located at distance I from the root. Let a = l/t. 
Then 



■log(m-i), _ 1 1 / {in-l)a^ 



/2it (r2 - a2)i/4 y q,2 + (m - 2f 



E 

.k=l 



+0 Q ,ifO<a<r. 
In addition, there is a constant c > 0, which depends only on m, such that 

We will prove this theorem in Section [3] 
Next, define 

p(a,i) =m(m-l)M-i-i|V^[„4](t)|2, (5) 



and 



t+Vt 

dt (6) 



p{a,T) = — p{a,t) 

The factor m{m — in ([5]) equals the number of vertices in the tree at the 

distance [at] > 1 from the root. Intuitively, p {a, T) is the average probability 
density of the event that we find a particle at the distance approximately at from 
the root if we measure its position at time approximately equal to T. Then, we have 
the following corollary of Theorem 1 1.51 

Corollary 1.6. Suppose that oji (a) ^ uj2 {a) . Then 

1 1 TflOp' 

lim p(a,r) = -===— -2, //■0<a<r, 

T^oo TT - oi^ j^{jn- 2) 

= 0, if a > r. 

The proof of this corollary is similar to the proof of Corollarv ll.2l and is omitted. 

The limit distribution is similar to the arcsine distribution (HJl, except it has a 
weighting factor, which is different for every m. A plot of the limit average distri- 
bution is shown in Figure |2] for m = 4 and m = 20. For the purposes of compar- 
ison we have additionally rescaled the support of the distribution so that supports 
are the same for both m. It can be seen that the walk on the tree of higher valency 
has higher density next to the border of the support. 

Note that the quantum walk on trees behaves (perhaps non-surprisingly) quite 
different from the classical random walk on trees. In the latter case, it is possible 
to show that the distribution of the particle distance from the root is asymptotically 
Gaussian with mean ct and standard deviation where c, a > (see |[T6l and 
ifTSl ). This result is very different from what we find in the quantum case. 
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Figure 2. The limit average probability distribution for the 
quantum walk on a homogeneous tree with valency m. The solid 
Une is for m = 4, and the dashed line is for m = 20. The support 
of the distribution is rescaled to [0, 1] interval. 




Figure 3. The transition amplitude of the return to the root for 
the nearest-neighbor quantum walk. The solid line is for integer 
lattice Z; the dash-dotted line is for 4-valent infinite tree T4. 

Consider now the amplitude of the return to zero at time t. 

Theorem 1.7. Let ipo (t) denote the transition amplitude of the return to the root at 
time tfor the nearest-neighbor quantum walk on the m-valent infinite tree started 
at the root. Suppose that m > 3. Then, for large t the following asymptotic ap- 
proximation is valid: 

^0 (t) = - Jt3/2"l2 [{2V^ t - n/i] + O (^) . 
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The plot of ipo (i) for m = Ais shown in Figure [3]by dash-dotted line. We can 
see that the frequency of oscillations in the return amplitude is higher than in the 
case m = 2. In addition, the absolute values of maxima decline faster. 

As a corollary of Theorem 11.71 we can see that the probability of the return to 
zero has the following asymptotic approximation: 

PO (t) = ""'^""'llT i [(2V^) t - vr/4] + O (t~y^) . 

IT [m — 2) ^ ^ 

If we compare these results with the case of the classical random waUc, we find 
two surprising facts. First, there is no exponential decay factor in the probability 
of return. The decay is polynomial of order t^^. Second the exponent in this poly- 
nomial decay does not depend on the valency of the tree although the frequency of 
oscillations and the overall constant does depend on it. 

The rest of the paper is organized as follows. Section |2] provides proofs for 
Theorems 11.11 and 11.41 concerning quantum walks on Z and Z'^. And Section [3] 
gives proofs for Theorems 11.51 and 11.71 concerning the nearest-neighbour quantum 
walk on homogeneous trees. 



2. Quantum walk on integer lattices 

Proof of Theorem ll.lt It is convenient to introduce additional notation. Let G 

be a graph {V, E) . For f,g^L'^ {V) , we define 

{f,g) = Y.J{x)9{x)- 
Then, for V = Z we can write 



fc=0 



01 kl 



where S^i is the shift operator that sends to (5fc_/. 

Let a[''^ := {X'')^^ . Then, it is easy to see that YlZ-oo 4''^ = (-2)^ )With 
(j){z) := Yli-L^i^^- Therefore, 



2m J z'+i 
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where the integration is over a small contour around 0. Hence, the transition am- 
plitude is given by the formula 



k 



jk) jit) ' 



k=0 



2m J z'+i 

f-27r 



A;! 

dz 



1 

"2^ 



e 



it{p{e)+ae)^Q 



where we made the change of variables z = e~*^ and where a = l/t. 

We can evaluate the asymptotic behavior of this integral by the method of sta- 
tionary phase (Chapter 4 in [TSl). The points of the stationary phase can be found 
from the equation: 

P' {6) = -a. (7) 

Suppose that this equations has K > solutions 9^ (a) . Then, the asymptotic 
contribution of the stationary point 9k to the integral above is given by the following 
formula: 

1 1 Jt{P{dk)+aek)±7ri/4 ^ Q I 



V27rt y^\P" {9k)\ \t. 
where the sign before ni/A depends on whether P" {9k) is positive or negative. By 
adding these contributions we obtain the first claim of the theorem. 

If the equation (O has no real solutions, then there are no points of stationary 
phase in interval [0, 27r) . In this case, we can apply the method of integration 
by parts. Usually, in this case the asymptotic approximation is of the order t~^. 
However, in our case it is smaller due to special properties of function P {9) and 
number a. 

Indeed, let us denote P (9) + a9 as fa {9) for shortness. Note that the first 
derivative is periodic with period 27r and therefore all other derivatives are also 
periodic with period 27r. In addition, if at is integer (which is exactly the case we 
consider, then exp {Ufa {9)) is periodic with period 27r. 

Since \ f' {9)\ / anywhere in interval [0, 27r) , hence we can use integration by 
parts in the following form: 

2n ^ i-2n ^ 



It Jo f {9) d9 V ; 

1 



1 ^gjt^(27r)\ _ 1 Cgit/c(0) 



/'(27r)V ; /'(O) 

git/a (e)^^ 



1 z-^- d 



it Jo d9\f'{9) 
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By using the special properties of the function /„ (9) , we can conclude that the 
first term is zero and therefore 

it Jo \f'{9)J 

In particular, this integral is O {t^^) ■ Since the function (1//' (9))' is periodic, 
hence the argument can be repeated . It is easy to see that it can be repeated 
indefinitely, and we obtain that the integral is less than c„t~" for every n. QED. 

Proof of Theorem ll.4t Let X be the adjacency matrix for 7? and H and V be 
the adjacency matrices that take into account only horizontal and vertical bonds, 
respectively. In other words, 

X = H + V, 

and 

Vij,kl = ^ik + ■ 

It is easy to see that 

That is, H and V commute. This implies that 

After we apply e**^ to ^ (0) = (5(0,0) ) the result at vertex (i, j) is Joi^j (*) where 
Sqi is the Kronecker delta and (t) is the wave function for the nearest-neighbor 
quantum walk at Z started at 0. Next, after we apply e**^, the result at vertex i,j 
is tpi (t) tjjj (t) since it equals the wave function of the nearest-neighbor walk on 
graph Zx (0, j) started with the initial data tpj (t) . QED. 

3. Nearest-neighbor quantum walk on trees 

Proof of Theorem ll.St 

First, note that 

fc=0 

k=0 

where (|w|) denotes the number of all possible paths with k edges that start at 
the root and end at vertex w. 
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Let Ak denote the number of paths from e to e that have length k and do not pass 
along a specific edge which is connected to e, say, do not pass along edge xi. Let 
Bk be the number of paths from e to e that have length k, without any additional 
restrictions. Let A{z) and B {z) denote the generating functions for Aj^ and B^, 
respectively, that is, 

oo oo 

A{z) = J2 Akz\ and B{z) = Y, Bkz\ 



k=0 



k=0 



where we set Aq = Bq = 1. 

From Lemma [3?T] proved below, it follows that 



J2ci+r{l)z' = A{zyB{z). 



r=0 



Hence, 



ci+r il) = 7r^<t> dz. 



27ri J z' 

where the integration is around a small circle around 0. 

Let \ w\ = I, then for the transition amplitude from e to w, we can write 



k=l 

oo 



k\ 



r=0 

which we can re- write as follows: 



(/ + r)! 2m 



i^i it) 



^jA{ziB{z)z^-^[Y: 



2iri 



{it/zf^ 

[i + ky. 



\k=0 



dz 



A{zyB{z) z' 



Jt/z 



l-l 
fc=0 



jit/z) 
k\ 



dz. 



The sum in the last Une gives zero contribution to the integral since neither A {z) 
nor B {z) has any singularity at 0. Hence, we can write 



i^i it) 



— i Aiz)^ B{z)z^-^e'*l'dz 
2m J 



1 

2TTi 

1 

2^i 



Ail/u) 



u 



'B{l/u) 



e^'^'du 
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where we used the substitution u = l/z, and 











u ' 












u 


U 



F{u) : 
G{u) : 

The second and third integrals are taken over a sufficiently large circle around the 
zero which includes all of the singularities of F (n) and G (n) . 

We calculate F (n) and G (u) explicitly below (Lemma [X2] and [33] ) . The func- 
tion G (n) is analytical at points u = ±m, therefore the only singularities of the 
integrand ai^e branch points of F{u) and G{u) at n = ±2y/m — 1. 

We want to find out the asymptotic approximation for those values of I which 
are comparable with t. Let I = at with a > 0. The we can write the transition 
amplitude as follows: 

tiji (t) = — /e^*["-^"'°g^(")lG(w)d'u. (8) 
27ri . 



Recall that r := — 1. Let us deform the contour of integration so that it goes 
first from — r to r just below the real axis, and then goes back just above the real 
axis. 

Let / (u) := u — ia log F{u). For real u G [— r, r] , we can compute Im/ (n) = 
{a/2) log (m — 1) , which is constant with respect to u. Hence, we can use the 
stationary phase approximation to this integral. 

In order to find the points of stationary phase, we need to solve the equation 
d (Re/ (u)) /du = 0. Since Im/ (n) is constant, it is the same as solving 

df (u) F' (u) 

du F [u] 

First, let us consider the case < a < r. 

For the part of the contour that lies in the upper part of the complex plane, we 
have: F (u) = (^u — i\/r'^ — u^^ / (2 (m — 1)) , hence /' (n) = 1 + a/\/r'^ — u'^ 

and equation ^ becomes —a = \Jr'^ — n?' which has no solutions in the inter- 
val (— r, r) for any positive a. Hence, the contribution of this part of the contour 
is asymptotically negligible provided that the integral along the other part of the 
contour has stationary points. 

For the part of the contour that Ues in the lower part of the complex plane, we 



have F [u) = + i-y/r^ — n^j / (2 (m — 1)) , and the equation ^ reduces to 
a = \/r'^ — V? ^ which has two solutions u\2 = i-y/r^ — o? for a <r. 
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Recall that the method of stationaiy phase says that if u is the only stationary 
point of function / (u) , located inside [a,b] , then 



J a 



2tt 



where the sign before iri/Ais positive if Re/" (n) > and negative if Re/" (n) < 
0. 



We compute F (tii,2) = {^^Vr"^ — o? + mj / (2 (m — 1)) . The second deriv- 
ative of / (u) can be evaluated at u\^2 as follows. 



— 



In addition, we have 



f («1,2) = T 2 



± (m — 2) \Jr'^ — o? 



ima 



and 



2(a2 + (m-2)2) 



/ m — 1 
a'^ + {m- 2y 



Hence, 

it) 



e 2 



log(m.-l)^ 



1 



1/2 



m — 1 



2^t (^2 _c,2)l/2 



a2 + (m - 2)' 



.k=i 



+ 0(1)). 



Here, the frequencies can be computed as 

a 



uji = a arctan 



+ \/r2 — 



— 



and 



u}2 = avr — uJi 
and the phases can be computed as 

ifi = — arctan 



ma 



(m - 2) Vr2 - a2 4 ' 



and 



y,2 = -TT - (yJi. 

Now, consider the case a > r. In this case, neither part of the contour has a point 
of stationary phase and for the large t, the boundary points of the interval [— r, r] 
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contiibute most to the integral. In this situation, we can estimate the integral by 
using integration by parts. Consider, for example, the integral 

hie) = f ' e''^^''^G{u)du, 

J~r+e 

where / (u) and G (u) are defined as continuous limits of the upper half-plane 
branches of / (u) and G (u) . Then we can write: 



hie) 



1 

it 



d 

du 



u) Giu) 



itf (u) 



r+s 



-r+£ 



fin) 
1 
it 



du 



Jtfiu) 



r+£ 



Giu) 
fin) 



du. (10) 



Since /' (u) = 1 + a/ \Jr'^ — u^, therefore — >• as it — )• itr, which implies 
that the first part of (fTOl) becomes zero as e — 0. For the second part, note that 

G' iu) = Aiu) + Biu)Vr^ - u^ + C (n) ^ 



for some functions A (u) , i? (n) , and C iu) analytic on [— r, r] , which implies that 

'G(m)V _ G' iu) Giu) „ 

I J 



f'iu)) fiu) [f'^u)] 



has singularities (u + r) and (r — u) at —r and r respectively. This im 



1/2 



plies that 



G{u) 
/'(«) 



is absolutely integrable at [— r, r] and therefore 



lim Ii ie) 

£-S>0 



~ t 



Giu) 



f'iu) 



du 



< -e 
~ t 



-(at/2) log(m-l) 



A similar estimate holds for the integral along the contour in the lower half- 
plane. This completes the proof of Theorem 1 1.5 1 

Here are the auxiliary results that we used in the proof. 

Lemma 3.1. Suppose that graph G is an infinite homogeneous tree with root e. 
Let Ak be the number of paths in G from etoe that have length k and do not pass 
along a specific edge which is connected to e. Let be the number of paths from e 
to e that have length k, without further restrictions, and let c^^ i\w\) be the number 
of paths from etow that have length k. Then, 



Ck [\w\ 



E 

ko+k\+...+k\, 
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Proof: Assume that each edge in the tree is oriented and has a label, x, which is 
chosen from the set {1, ... , m} . It is assumed that that the labels of edges around 
each vertex are all different. We write label x if we move in the direction of the 
orientation and x~^ if we move in the opposite direction. Let xixi^i...xi be the 
shortest path from e to w. There is a one-to-one correspondence between the set 
of shortest paths and vertices so we can write w = xixi-i...xi. Also, let Wi = 
XiXi-i...xi. This is one of the vertices on the shortest path from e to w. We write 
the edges in the path from right to left so that w\ is a neighbor of the root. 

Every path from e to w can be considered as the shortest path from e to w 
decorated with loops which can be attached at each of the points of the shortest 
path, Wi. In order to make sure that we do not double count the loops we forbid 
the loop attached at Wi to go along the edge that connects Wi to u^j+i. In this way, 
at every point of the path we know in which loop we are in: We are always in the 
loop attached at that ivi that has the largest length | Wi \ among all those vertices Wi 
that have akeady been visited. 

Let I = \w\ . The number of possible different loops that can be attached at wq, 
wi, wi-i is counted by Akg, ^fci, ■ ■ ■ , ^h^i, respectively, where k(),ki, . .., 
are the lengths of the loops. The number of different loops that can be attached 
atw = wi is counted by Bj^^ . Then, the total length of the path is ko + ki + + ki 
and by assumption it must be equal to k. Hence the total number of paths is 



Proof: The function B (z) is related to the Green function of the nearest-neighbor 
random walk on an infinite tree, which is well-known. (See Dynkin and Malyutov 
for the seminal contribution, and Lemma 1.24 on p. 9 in fVT\.) Hence, we can 
compute 



A-oA-i • • • A•,_l^fcr 



fco+fcl+...+fc^=fc 



QED. 



Lemma 3.2. 




(11) 



B{z) 



(m - 2) + my^l - 4 (m - 1) z'^ 
2(1- 7n2z2) 



It follows that 



G(z) 



(m — 2) z + m-sj z"^ — 4 {m — 1) 



(12) 



2 (z^ — 777,2^ 



QED 
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Note that we chose the branches of G (z) in such a way that the function is 
analytical outside the cut \—2\Jm — 1, 2\Jm — l] . In particular, this function does 
not have poles at itm. 

More precisely, the sign before the square root is determined by the rule that for 
sufficiently small t, 

. \/m — 1 



G{it) ^ -i- G C" 



m 



and 



^ / . N . \/m — 1 

G{-it) ^ I- G 

m 



Lemma 3.3. 

z \z J 2 (m — 1) 

Proof: In order to compute A (z) , we note that the following recursive relation 
holds. 

k-l 

A2k = (m - 1) ^ A2iA2{k-i-i)- (13) 

1=0 

Indeed, consider a path from e to e, that avoids the edge xi. There are m — 1 
possibilities to start the path. Suppose that the path starts with Xi, i / 1, so that 
the second point on the path is the endpoint of xi which we denote wi . Let r be 
the first time when the path returns to e. Then Wr-i = wi and the path from wi to 
Wr—i is one of the Ar-2 paths from wi to wi that avoid passing thi^ough the edge 
labelled Xj. The remainder of the path goes from e to e and it is one of the A2k~r 
paths that avoid the edge xi. The number r must be even, greater than and less 
than 2k. Hence we can write it as r = 2/ + 2, where < / < A; — 1. This implies 
the recursive formula (fT3] ). 

Next, we can use the recursion formula for Catalan numbers, 

k-l 

Ck = ^ ClCk-l-1, 
1=0 

and formula (fT3] ) in order to conclude that 

A2k = {m- l)^Cfc. 

By using the generating function for Catalan numbers, we obtain the following 
formula for ^ (z): 



It follows that 



i-yi 4 

^ ^ 2 (m - 1) ^2 

_ ^- -4(m- 1) 
2(m-l) 
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QED. 

The sign of the square root in the expression for F {z) is determined by the 
following rule: for all sufficiently small t, 



F {it) ^ -i/Vm - 1 G C", 

and 



F {-it) i/Vm - 1 G C+. 
Proof of Theorem ll.7t By ([8]l, we need to find asymptotics for 

^o{t) = ^.j>e''''G{u)du, (14) 

where 



— {m — 2)u + m-\A?" 
G{u) = — 



1 (tt^ — W?') 

and r = 2\Jm — 1. We can deform the contour so that it starts at — r, passes just 
below the real axis to r and then returns back to — r just above the real axis. Then, 
we find that 



1 r 



2 — /1, 2 



my — u 



4^0 {t) = - I 

The main contribution is produced by singular points itr. After integration by 
parts, we obtain the following formula. 



, , , 1 m 



1 -ue**" , r -2u^\ 



Vr'^ -vP'U^ - m"^ J-r {vP'-m?)'^ 



J.2 _ y^^ 



du+ / ^e""(in 



We can apply van der Corput's results (see JH, p. 24) to the first integral in the 
brackets and obtain the following asymptotic approximation. 



We can apply the integration by parts to the second integral and find that it is 
0{t^^) .It follows that 

VoW = ^,,„,T^^2 sin(rt-./4) + 0(t"^) 

1 ^^^-^^^'\:^^^2V^^^t-^/A)+0{t-^). 



0Ft3/2 _ 2)2 
QED. 
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